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METHOD FOR ESTIMATING PITCHING-MOMENT INTERFERENCE 

OF WING-BODY COMBIMTIONS AT SUPERSONIC SPEED 

By George E. Kaattari, Jack N. Nielsen, 
and William C. Pitts 


SUMMARY 


This report is the second in a series investigating the effects. of 
interference at supersonic speeds on the aerodynamic properties of wing- 
body combinations. The first report, MCA RM A51J04, 1951? presented a 
method of determining the lift of triangular, rectangular, and trape- 
zoidal wing-body combinations . The present report extends the same 
method to the determination of the pitching moments of these combinations. 
Calculated centers of pressure are compared with experimental values for 
nearly 100 triangular, rectangular, and trapezoidal wing-body combina- 
tions. The experimental and estimated locations of centers of pressure 
were found to correlate as follows: The correlation points for the tri- 

angular wing-body combinations have an average deviation of 0.008 body 
length about the line of best fit which is displaced 0.009 body length 
from the line of perfect agreement. The correlation points for the rec- 
tangular wing-body combinations have an average deviation of 0.015 body 
length about the line of best fit which is displaced 0.026 body length 
from the line of perfect agreement. The correlation points for the trape- 
zoidal wing-body combinations have an average deviation of 0.016 body 
length about the line of best fit which is displaced 0.017 body length 
from the line of perfect agreement. It is recommended for design pur- 
poses that these displacements of the lines of best fit from the line of 
perfect agreement be subtracted from the appropriate theoretical value 
for center-of -pres sure position. A numerical example illustrating the 
method is given. 


INTRODUCTION 


The first report of this series (reference l) reviewed the existing 
methods of calculating lift components of wing -body combinations . These 
methods are either laborious or restricted to particular wing-body 
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combinations . In order to fill the need for a simple, general method 
applicable to a variety of combinations, a generalization of the method 
of Nielsen, Kat^en, and Tang (reference 2) which successfully predicts 
the lift of triangular wing-body combinations is presented in refer- 
ence 1. The applicability of the method has been verified by comparison 
with experimental re suits involving a wide variety of supersonic wing- 
body combinations including nearly 100 triangular, rectangular, and trap- 
ezoidal wing-body combinations . The present report extends reference 1 
to .include a method for calculation of the centers of pressure of the 
lift components of a wing-body combination and thus permits the calcula- 
tion of the pitching moment of the combination. 
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SYMBOLS 

aspect ratio of wing panels joined together 
lift-curve slope based on exposed wing alone 

moment-curve slope based, on exposed wing -alone 

chord at wing-body juncture, inches 

-chord .at wing tip, inches • 

wing chord at spanwlse distance y from body axis, inches , _ 

body diameter, inches c 

ratio of lift of body in presence of wing exclusive of. nose to 

that of wing alone . " ' 1 . : : =: ' ‘ 

ratio of lift of body nose to that " of wing alone 

ratio of lift of wing in presence- of body to lift of wing alone 

lift force, pounds 

body length, inches . . “ ■ = : 

body length behind wing trailing edge, inches 

body length forward of wing leading edge, inches . . 

distance from most forward point, of body to center of pressure, t ‘ 

inches - 
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B(W) 


distance from most forward point of body to moment center of 
combination, inches 

distance from most forward point of body to center-of -pressure 
position of body nose, inches 

moment reference length, inches 

moment, pound-inches 

free-stream Mach number 

cotangent of leading-edge sweep angle 

free-stream dynamic pressure, pounds per square inch 

body radius, inches 

exposed wing area, square inches 

semispan of wing-body combination, inches 

volume of body considering the body as cylindrical behind the 
position of maximum cross section, cubic inches 

streamwise, spanwise, and vertical coordinates, respectively 

chordwise distance from leading edge of wing-body juncture to 
center of pressure, inches 

chordwise distance from leading edge of wing-body juncture to 
local center of pressure at spanwise location y, inches 

angle of attack of body 

local angle of attack at spanwise distance y from body axis 
VMq 2 - 1 

leading- edge sweep angle 
taper ratio 



Subscripts 


body alone 

body in presence of wing minus body nose 
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C wing -body combination 

C-N combination minus nose 

N nose of combination 

W wing alone 

W(B) wing in presence of body 


Superscripts 

s slender body 

u upwasb theory 


ANALYSIS 


Since slender -body results are used in this analysis, the method is 
restricted to those wing-body combinations for which slender-body results 
are available. This means that at the present time swept-forward lead- 
ing edges or swept -back trailing edges are generally precluded. A typi- 
cal configuration is shown in figure 1, which gives the necessary 
dimensions required in lift and moment calculations. 

The moment of a combination is the summation of the moments of all 
lift forces acting on the combination. To compute these moments requires 
that all lift forces and their locations (centers of pressure) be found. 
The methods of determining these forces and of finding their correspond- 
ing locations are now discussed. 


Lift Components 


The lift components of a combination are considered to be the lift 
of the body nose 1 L^, the lift of the wing in presence of the body 
and the additional lift of the body due to the presence of the wing 
(exclusive of the nose) These components, following reference 1, 

1 The nose of the body is that part of the body in front of the wing 
leading-edge body juncture, or, in cases where the wing is mounted on 
an expanding body section, is that part of the body in' front of the 
position of maximum cross section. • 
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are conveniently evaluated by multiplying the lift of the wing alone 
(defined as the exposed wing panels joined together) by appropriate fac- 
tors %, %( B ), and K B ( W ). Thus with L N = %L W , Iy( B ) = %( B )L, 
and I>B(y) = the lift °f the combination is given by 

^ = E% + %( B ) + %(w)3 L W (!) 

where . may be either the theoretical or experimental value. The 
value for calculated by linear theory is used throughout this 

report . 

Evaluation of %.- The factor Ktf is defined as 



( 2 ) 


giving the ratio of nose lift to that of the wing alone. The lift of 
the nose in equation (2) may be evaluated by slender-body theory 


l N „ - 

— = 2rcr 2 
qa. 


(3) 


The slender-body value for Ljj was used in the present report and in 
reference 1. It is natural, however, to expect improved accuracy by 
using linear theory, and linear- theory results for Ljj may be used if 
desired. 


Evaluation of 




body, the factor %( B ) is 


For lift of the wing in the presence of the 
defined as follows : 


%(B) “ 


^(B) 

% 


00 


It was shown in reference 1 that the ratio of the lift of a wing in com- 
bination with an infinite cyli n d r ical body to that of wing alone as given 

by slender-body theory ^j(b)^ S ^ is "vslfd even for high aspect-ratio 

wings. The value of is presented in figure 2 as a function 

°f r / s m 8Jid has been used for ^(b) i Q this report. 


Evaluation of ^b(w) *“ 'fh e method of determining the additional 

lift of the body due to the presence of the wing, as given in reference 1, 
is based on the concept of the body supporting lifting disturbances 
emanating from the wing. The region between, the Mach helices on the 
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body originating at the wing leading and trailing edges, as shown in 
figure 3( a )> represents the region supporting most of the lift .carried 
over to. the body from the wing. Simplification of this nonpla n a r model 
to the equivalent planar representation shown in figure 3(b) is desirable 
for purposes of calculation. The body is then represented by a flat 
plate at zero angle of attack (ccg = 0) and the Mach helices of figure 3( a ) 
become Mach lines in figure 3(b). 


The value of lift carried over onto the body by a half wing with 
supersonic leading edge is given by 


L = 


gm 


=/V 

■1 0-0 Of, 


C r +T1 


cos 


.-1 p 


I- + pmn 


TJ + m | 


d£ 


(5) 


with the coordinate system of figure 3(b). (See reference 1.) This 
result, doubled to account for the lift of two half wings and divided by 
the lift of the wing alone, gives %(w) * • . 


K B(W) 


8gm . 

.c/eV-i (^) (lf)(v - 1 ) ( pCl «) w 


t 

(gm+l) |^ + gm 

s 

cos” 1 

1 + (l+gm)g A 
c r 

\ gm+V 

gm 

gm + (gm+l) ^ 






«/g 2 m 2 -l 

( gm+iy 



<s/ g 2 m 2 -l 
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( 6 ) 
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when mg>l. 
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Similarly for subsonic edges there is obtained 

: />f 


L = 


SqayCpm) 3 / 2 P^ p c r +T l J $ 

jrp(pm+l) J 0 


*/ ml+Tj 


dl 


giving 
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p 2 !! 2 


(l+Pm)- 


* (1+A - ) (i) (v - 1 ) (p°i«) 


w 


pm + (l+pm)|^ 


Pm 


3/2 


Pm + (l+pm)|^ 
‘-r 

1/2 

(l+pm)M 

2 

_ Q _ 

tanh -1 / ^ 

pin. 

“ C “ 

pm 

— — 

V Pm + (-l+pm)|| 


(T) 


(8) 


when mp< 1. 


Tip effects are not considered and the analysis is confined to cases in 
which the Mach line from the leading edge of the wing tip falls behind 
the assumed region of lift carried over onto the body. This condition - 
Imposes the restriction 


PA (l+X) Qj + l) > k (9) 

on the wings for which equations (6) and (8) are to apply. Equations (6) 
and (8) are represented graphically by plotting the quantity 

%(W) ( 1+ *> (v - 1 ) (s°La) 

W 

as a function of pd/er for constant values of mp in figure 4, which 
is to serve as a design chart in determining %(w) subject to the 
restriction of equation (9). 


In reference 1, the following selection rule was given for choosing 
between %(y) and If PA(l+l) +1^ <k, use the slender- 


body value of from, figure 2. If pA(l+\) +1^—4, use the 
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value of K B ( W ) given by figure^ 4 unless k b(w) > k b(w) < ' S ? • In such 
cases use 


Center of Pressure Positions 


Center- of pressure of the body nose .- The location of center of 
pressure of the body nose as given by slender-body theory is 



( 10 ) 


where V is the volume of the body considering the body as cylindrical 
behind the position of maximum cross section. Equation (10) is used, in 
this report, although the linear-theory result for 2 N , if available, 
may be expected to give improved accuracy. 

Center of pressure of wing in presence of body .- Several methods of 
determining the center of pressure of a wing in the presence of a body 
will now be presented. ; 

The, center of pressure of a triangular wing in the presence of an 
infinite cylindrical body as given by slender-body theory (reference 3), 
in percent of the exposed wing root chord measured from the leading edge 
of the wing-body juncture, is 



An alternate method for evaluating center -of -pres sure location of a 
triangular wing-body combination is to suppose that the exposed wings 
are operating in the upwash field of the body alone and then to calcu- 
late the resultant center-of -pressure location using strip theory. Neg- 
lecting any effect of the nose, it has been pointed out (reference 4) 
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that the upflow angle due to the body varies spanwise on the horizontal 
plane of symmetry as 


Oy = Ct£ 



( 12 ) 


where y is the lateral distance from the body axis. The wing is thus 
effectively twisted by the body -alone flow. If now the upwash angle 
given by equation (12) is taken into account by using strip theory, an 
approximate value of lift is given as 


L W(B) (u) « “ °y c y d ; 


P ^ r 

The moment about the leading edge of the root chord is 

' s m 


(13) 


/ s 1 r B m _ 

%(b)^ uj =-i 7 “y x y c y d y 




on the assumption that the center of pressure of the strip is at the 
midchord. Dividing moment by lift then gives for the center -of -pressure 
location for a triangular wing-body combination 




(t)' 


2 + 


(15) 


The results of equations (ll) and (15) are presented in figure 5 as a 
function of r/sm. In addition, the value of center of pressure of the 
wing alone as determined by linear theory is indicated. It is signifi- 
cant that all three methods give essentially the same result for the 
center-of -pressure location of the wing in presence of the body. It may 
be concluded that (x/ c r )^ for wing alone (defined as exposed wing 
panels joined together), although independent of r/s m , gives a suffi- 
ciently accurate representation of (x/c r for triangular wings in 
presence of the body. 

If slender-body theory is applied to rectangular wings in combina- 
tion, the erroneous result is obtained that all lift, and therefore the 
center of pressure, is at the wing leading edge. While this result is 
valid for vanishing aspect ratio, it is obviously not valid in general. 
On the other hand, by strip theory, the center of pressure is given at 
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the midchord and 'is independent of the aspect ratio. This value is 
exact only in the case of. vanishing chord and is approximately true for 
moderate to high aspect ratios. The center-of -pressure location of wing 
alone as predicted by linear theory exhibits a shift toward the leading 
edge from the midchord position with decreasing aspect ratio. 


x_) _ 3PA-2 

c r^‘ 6pA-3 * * * 


(16) 


Equation (l6) is valid for pA>l. For f3A<l, negative lifting pres- 
sures due to tip effects develop on rearward areas of the wing, moving 
the center of pressure nearer the wing leading edge. Thus, wing-alone 
center-of -pressure location as predicted by linear theory approaches the 
value given by strip theory for wings (in presence of body) of high 
aspect ratios and shows a location more in accordance with slender-body- 
theory results at low aspect ratios. It. is therefore concluded that the 
center of pressure of the wing alone for all aspect ratios is more repre- 
sentative of the center of pressure of the lift on the wing in presence 
of a body than the result given by either slender-body theory or strip 
theory. 


For trapezoidal wings of no trailing-edge Bweep, slender-body theory 
gives all the lift, and hence center of pressure, on the portions of the 
wing forward of the leading edge of the tip chord. In general, h owe vet, 
lift is known to exist over the entire wing and the slender-body result 
for center-of -pressure location is too far forward at high aspect ratios. 
Strip theory, on the other hand, principally by not accounting for tip 
effects, generally gives a center -of -pressure location too far -aft of 
the wing leading edge particularly at low aspect ratios. For large 
aspect ratios wing-alone theory is in good accord with strip theory, and 
at low aspect ratios, with slender-body theory. Since strip theory is 
reliable only at high aspect ratios, it can be concluded that wing-alone' 
theory is best for the entire aspect-ratio range. 

On the basis of the foregoing comparison of wing-alone theory with 
slender-body theory and strip theory for triangular, rectangular, and ' 
trapezoidal wings in combination with a body, it is concluded that of 
these three theories wing-alone theory is the best for representing the 
center of pressure of the exposed wing panels thro ugh out the aspect-ratio 
range; Some simple charts to assist in estimating these center-of- 
pressure positions are now presented. 

Figures 6(a), 6(b), and 6(c) give the variation of (x/c r )^ with 
pA for wings of no trailing-edge sweep, no midchord sweep, and no 
leading-edge sweep, respectively, for taper. ratios of X = 0, l/2, 
and 1. Tbe curves giving (x/c r )^ are extrapolated to the l imi ting 
values given by slender-body theory at |3A =-0, for Which case 
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slender-body theory is valid. Figures 6(a), 6(b), and 6(c) serve as 
design charts, and the value of (x/c r )y(£) for any given wing of this 
family may be found by suitable interpolation. If more accurate results 
are available for the wing-alone center of pressure, these results may 
be used instead of those found from figure 6. 


Center of pressure of lift on body due to the wing .- The method of 
determining center-of -pressure location of the lift on the body due to 
the wing is based on the same model used to determine the lift in refer- 
ence 1. (See fig. 3*) The moment of the lift (equation (5)) carried 
onto the body by a wing with a supersonic leading edge is 


*%(¥) " 




p 2 *-/? 


r(3m n 

== / ' an / 

m 2 -l^o ^ti 


£ 


c r +Ti - 5 - + PmTJ 

£ cos" 1 


tj + m| 


d| 


(IT) 


with the coordinate system of figure 3(h). This result, doubled to 
account for the lift of two half -wings, gives 


^(W) = 


3* 


c r S ( 


F 


2pd 


2mp+5 ( pa/cr (pd/c r ) 2 1 


r L 3(mp+l) 3 (mp+1) 


*/m 2 p £ 


-.3 


at (Pd/cr) 3 

■ 1 

r / m 2 p a 

(1+mp) 2 _ 

(i+Sf) ■ 

2mP+5 

V pay 

L 3(mP+l) S J 


COS 


-1 


Pm 


1 + M (mp+1) 


mp + M (mp+1) 

Cr 


l- (n_T 1 

Vmp+l J i 

J - cos 1 ~ 


Vm 2 p 2 -] 


mp 


(18) 


The center -of -pres sure location is then found using 
tion (6) and the moment from equation (l8) as follows: 


from equa- 


G 


A') 


Cr/ B(W) 


%(W) 

^(w)^ 


%(W) 


K B(w) L W <i r 


(19) 
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Similarly for wings with subsonic edges there is obtained 

II 

' ,c r. +T l p 


^(W) " 


QqayCmg) 3/f£ rd- ■ n£ r +T l | J 3 " ^ 

113 s (pm+l) *-'0 v -'q */m£+T] 


d| 


giving 


( 20 ) 


4qa« 

%(W) = — °r 3 < 


/m 2 3 2 +m3 (m3+l 
v c r 


9m3(m3+l)'' 


(8m3+2^-)m 2 3 s +(l^mP+6)m3^ : + 


3(mP-3)(m P+ X) 2 (M f ]- • < a “P+ 21t >” 3 f . (^losh- 1 /— — jf 

\ c r/ J 9m3(m3+l) s 3m|3 W/ J (mp+1) £d 


X J 
( 21 ) 


The moment of equation (21) with Kg( w \ of equation (8) is used in 
equation (19) to give the center of pressure of the lift Qn an infinite 
cylindrical body due to the wing. The results for center of pressure 
for both supersonic and subsonic cases are presented as a function 
of 3&/cr with m3 as the parameter in figure 7> which is to serve as 
a design chart. It is notable that the effect of m3 is small. 


According to the results of slender-body theory, the location of 
the center of pressure of the body in presence of the wing measured aft 
of the leading edge - of the wing-body juncture is 



where Kg(y}( 8 ) and K^ B )( S ) are functions of r/sm (fig. 2) 

/ x \( B ) 

( cj - J, is give 11 by equation (ll). Equation (22) is also plotted 
^ 'w(B) 


and 
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against pd/e r with mg as the parameter in figure 7 for comparison 
with the previous analysis. Again, m(3 is seen to have no large effect 

/_ \(b) 

on the center of pressure. The slender-body value of ( — ) is seen 

^^(w) 

to begin at a value of 0.5 for pd/e r =0 in agreement with the results 
of the previous analysis. However, for increasing values of gd/c r the 

/- \(b) p 

slender-body results show an asymptotic TnaviTnirm of ( — ) 

v ° r 4(w) 

by the previous analysis, which accounts for afterbody, the value 


= — : whereas 
3 


of 




' ry B( W) 


continually increases with M;. The latter result is to 


be anticipated by a consideration of figure 3(a)* For a given geometry, 
an increase in Mach number causes a primary portion of the pressure dis- 
turbance carried onto the body to sweep beyond the wing trailing edge. 
Similarly, a decrease in chord with a given Mach number and body diameter 
moves the wing trailing edge ahead of the primary portion of the lift 
disturbance carried onto the body. In view of the foregoing arguments 
and since slender-body theory does not properly account for an afterbody. 


the present method of determining 

with afterbodies. Figure 7 serves 
combinations with afterbodies. 


(f) 


is applied to combinations 

'B(W) 

as a design chart for all wing-body 


Center of pressure of the combination .- The center-of -pressure loca- 
tion of a complete configuration in terms of body length is given by 


fl \ = %% + %(b) z w(b) + %(w)%(w) 
1 c *[% + %(b) + %(¥)] 


(23) 


where Zjj is given by equation (10) and 


Z W(B) “ l f + (Jr) Cr 

r W(B) 

(2V) 

^B(W) = l f + \JZ) c r 

X r/ B(W) 

(25) 


V. 
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NUMERICAL EXAMPLE 


To illustrate the use of the method developed in the foregoing sec- 
tions, the determination of the center-of -pressure location of a trape- 
zoidal wing-body combination is now presented. , Given that c^ = 1.500, 
c r = 3.878, r = 0.850, Sin = 3.790, Mq = 2.87, Z f = 16.66, l = 24.00, 

V = 39*98 and no midchord sweep, the basic parameters may be evaluated 
as follows : 


A = 


4(2.940) 
1.5 + 3.878 


= 2.19, aspect ratio of the wing alone 


p = < i/m 0 2 -1 = J 2.87 a -l = 2.69 


PA = 5.89^ effective aspect ratio 


r/s m 

X 

m 


0.224 body radius, semispan ratio 

- = O.387 taper ratio 
3.878 

2.47 

3.878 - 1.5 


m|3 = 6.64 


M = 1.18 

c r 3.878 


The value of (3 




from the charts of reference 5 is 3.85* 


The value of the parameter "in equation (9) is 


e A < 1+X > Qf + X ) - (5-89K1.387) + l) >b 


The preceding information now enables the lift parameters in equa- 
tion (23) to be determined as follows: 

K n given by equations (2) and (3) is 


% " 


2jtr‘ 


2fl(Q.85)‘ 


(%.-*)(<*«*) (VO /p 2.9M1.50 + 3.878)( 3 .85)/2.69 


= 0.200 
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%(B)( S ) from figure 3, for r/s m = 0.224, is 

%(B) = 1,18 

■%(W) f rom figure 4 in parametric form is 


k B(W) 


(poie) -a + x) (f! - l) = fcAl 


The value of is thus 


4.4l 


K B ( W ) = 0.24 

V ; (3*85) ( 1 * 387 ) ( 3 * 46 ) 

Now determine K B ( W ) ( B ) from figure 2: 

k b(w)^ = °*31 

Since %(^) < K b(w) ^ ’ fciie 781116 of %(-w) is to be used. 

The remaining parameters 7 N , 7 W ( B ), and 7 B ( ¥ ) are obtained in the 
following manner: ' ' - 1 - 1 . 

from equation ( 10 ) is 

T„ .1(2.. -J L.) - 2k [ 1 - , 39-96 1 , 6 .h 

15 V cr 2 zy L it(a 85 ) 2 ( 2 t) J 

1 W ( B ) from equation (2b) is 

h(B) = *f + {£) °r = 16.06 + (0.^9) ( 3 . 787 ) = 17.96 

_ Mb) 

where the value of f J = 0.49 was obtained by interpolation from 
figure 6(b). ^ (B) 

from equation ( 25 ) is 

*B(W) = Z f + ( #- ) c r = 16.06 + (1.08) (3.878) = 20.25 
v r y B(W) 


where the value for 


( — ) , from figure 7, is 1 . 08 . 



1 6 


NACA KM A52B06 


f 



The distance to the center of pressure of the 'combination from the 
most forward point of the body is then given by 


M = + Ky(B)?W(B) + Kb(W)*B(W) 

\2 J Q l[% + %(b) + %(w)3 


(23) 


(0.200)(6.4) + (1.18)(17.96) + (0.24) (20.25) 
2k (0.200 + 1.18 + 0.2k ) 


0.702 


EXPERIMENTAL VERIFICATION 


The foregoing analysis has been applied to the calculation of the 
centers of ■ pressure of nearly 100 wing-body combinations of widely vary- 
ing plan form. The results are compared with the experimental center of 
pressure found by putting the experimental values of and Cj/^ into 

the expression 



where Z r is the moment reference length in inches. Tables I, II, 
and III summarize the experimental and calculated centers of pressure 
for the triangular, rectangular, and trapezoidal wing-body combinations 
and also give the theoretical centers of pressure for the lift of the 
body nose, the lift of the wing in presence of the body, and 'of the 
additional lift on the body due to the presence of the wing. The corre- 
lation between the experimental and estimated results for the triangular, 
rectangular, and trapezoidal wing-body combinations are shown in fig- 
ures 8, 9, and 10, respectively. 

The sources of the test data are listed in references 6 to 27; some 
of the test data are unpublished. Seme difficulty was met in trying to 
determine values of lift- and moment -curve slopes from published curves 
because of irregularities in the data, in such cases the average values 
of and over a small angle-of -attack 'range about a = 0 were 

used. Furthermore, some experimental lift and moment-curve slopes were 
questionable. In one case, data on similar configurations from differ- 
ent facilities (of different Reynolds numbers') gave a difference of the 
order of 10 percent - in the lift-curve slopes. In view of the foregoing, 
some of the deviation found in the correlations of experimental and cal- 
culated center -of -pres sure positions can be ascribed to questionable 
experimental data . 
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Since the division of lift and moment between wing and body was 
not given by experiment , comparison with the method for center of pres- 
sure for the lift of the body nose, the lift of the wing in presence of 
body, and the additional lift on the body in presence of the wing could 
not be made directly. The correlation was therefore made on the basis 
of center of pressure of the entire combination. 


Triangular Wing-Body Combinations 


The correlation of est ima ted center-of -pressure positions with those 
of experiment for triangular wing-body combinations is presented in fig- 
ure 8. Included in figure 8 is a line of perfect agreement and also 
lines of ±0.05 1 deviation from perfect agreement. The" correlation of 
estimated values with those of experiment is excellent; the line of mean 
correlation is displaced 0.009 body length from the line of perfect agree- 
ment and the average deviation from, the mean is only 0.008 body length. 


Rectangular Wing-Body Combinations 


The correlation of estimated center -of -pres sure positions with those 
of experiment for rectangular wing-body combinations Is presented in 
figure 9* The correlation of estimated values with those of experiment 
is not so good as that for triangular wing-body combinations. The line 
of mean correlation is displaced about 0.026 body length from the line 
of perfect agreement. The average deviation from the mean is 0.015 body 
length. A possible effect that can explain the difference in correla- 
tion between the triangular and rectangular wing-body combinations is 
that due to the wing tip. It may be seen that the lift carry-over from 
a rectangular wing onto the shaded area of the body shown in figure 3(b) 
is independent of span, provided that (3A ^2, and may be considered that 
due to an infinite wing. In order to form a finite wing, a "canceling 
wing" must be superposed on the infinite wing to form a wing tip. This 
canceling wing generates a negative lift which carries over in part onto 
the body aft of the trailing edge of the wing at a distance which depends 
primarily on the Mach number and wing semispan. While this negative lift 
carry-over is probably small, its effect on the over-all moment and 
center-of -pressure position of the combination may be appreciable due to 
the large moment arm involved. Since no account was taken of this 
decreased lift on the afterbody, the calculated centers of pressure for 
the rectangular wing-body combinations are too far aft. Triangular wings, 
having no tip chord, may be expected to have less wing-tip effects than 
rectangular wings. 
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Trapezoidal Wing-Body Combinations 


The correlation of estimated center-of -pressure positions with those 
of experiment for trapezoidal wing-body combinations is presented in 
figure 10. Correlation of estimated values with those of experiment is 
generally within the ±0.05 1 correlation "limits excluding the combina- 
tions with no afterbody. The line of mean correlation is displaced 
0.017 body length from perfect agreement and the average deviation from 
the mean is 0.016 body length. It is notable. that the mean displacement 
from correlation for trapezoidal wing-body combinations is intermediate 
between those of triangular and rectangular wing-body combinations as 
might be surmised. 


Combinations With No Afterbody 


The present equations are applicable only to those combinations 
having an afterbody. .However, as a matter of interest, some combinations 
with no afterbody, the centers of pressure of which were computed on the 
basis of having an afterbody, are presented in figures 8, 9, and 10 . 
(denoted by flagged symbols ) . A sufficient number is not included for 
each type of combination to make a reliable statistical comparison. How- 
ever, on the whole these combinations show a somewhat greater mean dis- 
placement of the line of best fit from the line of perfect agreement 
than do combinations with an afterbody, indicating that an afterbody has 
an appreciable effect. 


DISCUSSION 


The question arises how large an error can be made by neglecting 
interference in determining the center of pressure of a combination. 

The answer depends on how the components of the combination are added 
together neglecting interference. First, the combination may be assumed . 
to act as the sum of the body alone plus the exposed wings joined 
together. On the other hand, the combination may be 'considered as the 
Siam of the body alone and the entire wing where the entire wing is formed 
by the extension of the wing leading and trailing edges to the body cen- 
ter line. For the. rectangular wing-body combinations considered, estimat- 
ing the center-of -pressure positions using either approach gives as good 
a correlation as is given by the present method. For triangular wing- 
body combinations the center-of -pres sure locations computed on the basis 
of the exposed wing, neglecting interference, do not correlate as well, 
as those computed by the" present method. The correlation scatters from 
near agreement up to values of about 20 percent of the body length 
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forward of the experimental values. The scatter of the correlation 
points becomes even greater when the center of pressure is computed on 
the basis of the entire wing. For trapezoidal wing-body combinations, 
a better correlation is attained than in the case of triangular wings; 
however, the correlation by either approach is not so good as with the - 
present method showing a greater scatter about the line of best fit, and, 
in general, giving eenter-of -pressure positions too far forward on the 
combinations. 

It is recommended that the present method be used for calculating 
the pitching moment of a combination for the following reasons: First, 

the present method is rational and is applicable to a variety of wing- 
body combinations. Second, both the center-of-pressure position and the 
lift are provided in determining the pitching moment of a combination. 
Third, the breakdown of the combination moment into its components iB 
given. 

For purposes of calculation, it is recommended that the appropriate 
displacement value of the lin e of best fit from the line of perfect 
agreement be applied as a correction in determining the center-of- 
pressure position of a combination by subtracting it from the center-of- 
pressure position given by theory. When this correction Is applied, an 
average deviation of 0.016 body length or less, depending on the wing 
plan form, may be expected In .the location of center-of-pressure posi- 
tion. The uncertainty in the center-of-pressure position in terms of 
the wing mean aerodynamic chord Is, of course, larger and depends on the 
ratio of body length to mean aerodynamic chord. 


CONCLUSIONS 


On the basis of the correlations between the estimated and experi- 
mental center-of-pressure locations presented. for nearly 100 triangular, 
rectangular, and trapezoidal wing-body; Mach wave configurations. It was 
found that by the method of this report the following limits of correla- 
tion were attained: The correlation points for the triangular wing-body 

combinations have an average deviation of 0.008 body length about the 
line of best fit which is displaced 0.009 body length from the line of 
perfect agreement. The correlation points for the rectangular wing -body 
combinations have an average deviation of 0.015 body length about the 
line of best fit which, is displaced 0.026 body length from the line of 
perfect agreement. The correlation points for the trapezoidal wing-body 
combinations have an average deviation of 0.0l6 body length about the 
line of best fit which is displaced 0.017 body length from the line of 
perfect agreement. It 1 b recommended for design purposes that the dis- 
placement of the line of best fit from the line of perfect agreement be 
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subtracted from the appropriate theoretical value for center-of -pressure 
position. ' 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. 
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TABLE I.- SUMMARY OF AERODYNAMIC AND GEOMETRIC CHARACTERISTICS 
AND TEST CONDITIONS FOR TRIANGULAR WING BODY COMBINATIONS 
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*09 

*25 

1.18 

2.51 

.23 

3*83 

n-Bo 

S3.k6 

26. kk 

27.W 

.k8 

7-S3 

— 

3-95 

*kT 

7*k5 

26 

6-3 

-10 

•25 

IJfl 

2.83 

*27 

k.31 

11.00 

28.k6 

28.72 

27.k8 

•k8 

T.63 

— 

k.kl 

.k6 

T-ko 

27 

6.7 

.10 

.25 

l-ifl 

3*17 

.32 

k.8k 

11.80 

28.k6 

29.01 

27.k5 

.k8 

7.63 

— 

k.88 

.k6 

7-39 

28 

7-2 

.11 

.2k 

l.lB 

3.53 

• 38 

5-38 

11.83 

26. k6 

29. ko 

27.kk 

.kd 

7.6e 

— 

5.36 

•k6 

7-35 

29 

6.7 

.15 

*29 

1.21 

3*10 

*k8 

5*13 

11.B0 

55.29 

55.79 

51.31 

.90 

ik.25 

— 

k-91 

•90 

lk.26 

30 

7*9 

.18 

-27 

1.21 

3.96 

.71 

8.58 

u.8o 

55.89 

58.29 

50.76 

.89 

ik.09 

— 

6.21 

•87 

13-85 

31 

9.5 

.25 

•29 

1.21 

k.00 

1.00 

7-00 

11.80 

55.29 

58.87 

kS.Sk 

.86 

13.70 

— 

7.26 

.85 

13.50 

32 

7-8 

.58 


L.33 

3-92 

2.29 

913 

1-95 

k.15 

k.59 

3.70 

Jte 

5*29 

2.k8 

8.93 

.kl 

5-00 

33 

6.9 

•52 


1.33 

3-39 

1-77 

7*77 

1*95 

kug 

k.W 

3.73 

-ka 

5.3k 

1.72 

7.68 

.kL 

5.16 
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TABLE II.- SUMMARY OF AERODYNAMIC AND GEOMETRIC CHARACTERISTICS 
AND TEST COEDITIONS FOR RECTANGULAR WING-BODY COMBINATIONS 


Ho. 

Sketch 

"o 

R 

Wing . 
Section 


PA 

r 

Sm 

1 

B 

Refer- 

ence 

Facility 

1 

■ 

1.93 

o.l9xio s 

l b.o. 

0.59 

9.41 

0.172 

12.4 

n.o 

24 

Iangley 
9 in. 

2 

j— 

1.93 

,19x10 s 


.59 

7.35 

.210 

12.4 

n.o 

24 

Langley 
9 in. 

3 

4" 

1.93' 

.19x10 s 

■SB 

.59 

5.24 

.273 

12.4 

11.0 

24 

Langley 
9 in. 

4 


1.93 

,19x10 s 


• 59 

3.16" 

.382 

12.4 

11.0 

24 

mm 

5 

-f— 

1.92 

.4oxlo 6 

gg 

1.25 

5.64 

.140 

11.2 

10.2 

25 

Langley 
9 in. 

6 

'»■ i— 

1.28 

. 56 xlo e 

®hex. 

1.30 

4.27 

.153 

8.8 

13.1 

19 

Aberdeen 

7 

' i 

1.28 

1.12x10 s 

hex. 

2.60 

B2| 

.265 


12.2 

19 

Aberdeen 

8 

■■'I- 

1.72 

1.12x10 s 

hex. 

2.60 

1.87 

.265 


12.2 

20 

Aberdeen 

9 

■■■■I— 

1.72 

. 56 x 10 ® 

hex. 

1.30 

- 

7.48 

.153 

8.8 

13.1 

20 

Aberdeen 

10 

H 

1.72 

.66x10® 

hex. 

1.54 

1.87 



21.5 

20 

Aberdeen 

11 


1.62 

.40x10® 

— 

1.25 

1.66 

.350 

6.8 

9.8 


Iangley 
9 in. 

12 

1 II 

1.93 

,40x10 s 

- - - 

1.25 

2.14 

.350 

6.8 

9.3 


Langley 
9 in. 

13 


2.40 

.40x10® 

— 

1.25 

2,84 

.350 

6.8 

9.8 


Langley 
9 in. 

14 


1.90 

. 50 x 10 s 

3 d.v. 

1.47 

1.86 

.384 

0.0 

14.5 

27 

Mich. V. 

15 

-| 

2.00 

.79x10® 

d.w. 

3.00 

4.76 

cn 

CO 

0 


10.9 


Ames 
1x3 ft 

16 . 

-1 

1.50 

. 91 x 10 ® 

1 

ESI 

3.08 

.083 


10.9 


Ames 
1x3 ft 

17 

■ | ■ 

1.93 

.18x10? 

| - 

'.59 

3.17. 

. .382 

12.4 

11.0 


Langley 
9 in. 

18 

■ 

1.62 

.21x10® 


.59 

2 . 45 . 

CVJ 

CO 

m 

12.4 

n.o 


Tangley 


9 in." 

19. 


2.00 

— 

— 

1.32 

1.73 

.333 


24 ;0 

26 

Dainger— 

field 

,20 

■■ 'f 

2.00 

_ _ _ 

. 

1.32 

3.46 

.200 


24.0 

26 

Dainger— 

field 


^b.o. Indicates biconyex 
%ex. indicates hexagonal 
3 d.v. indicates double wedge 
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TABLE II.- CONCLUDED 


No. 

• 


B 

% 


Theoretical 

Expert ■enteJL [ 

Lift 


Lift 

cent* 

ES 

ar of [ 
issure [ 

tjsgm 

BSSfl 

0259 

mm 

EC1 


■a 


[jgg 

B251 

msSi 



1 

lfl.8 

0.17 

0.12 

UA 

3-79 

0.6K 

5.W 

1.68 

*.12 

V.65 

3.87 

o.w 

3.34 

o.6v 

3.44 

O.H-l 

5.10 

2 

14.7 

.22 

.16 

1.17 

3.73 

.82 

5.76 

1.68 

V.12 

4.63 

3.83 

.vv 

5*47 

.62 

5-47 

.Vl 

3-U 

3 

10.5 

.32 

.22 

1.23 

3.62 

1-15 

6. VI 

1.68 

V.12 

4.65 

3.74 

.43 

5.33 

1.15 

6.69 

•42 

5-25 

V 

6.3 

.5 6 

.*0 

1.33 

3.3T 

1.90 

T.73 

1.68 

V.ll 

4.S3 

3.61 

.Vl 

5.15 

1.90 

7.16 

-36 

4.36 

5 

U-3 

.06 

.12 

1.11 

3.65 

-23 

V.71 

1.68 

V.12 

4.78 

v.u 

•47 

3.87 

•23 

V.37 

.46 

5.75 

6 

a.j 

.06 

.14 

1.12 

3.53 

.22 

4.67. 

3.18 

8.80 

9.37 

8.61 

.57 

6.89 

.19 

5.06 

.36 

6.78 

7 

2.1 

.10 

.38 

1.22 

2.13 

.22 

3.63 

3.18 

8.56 

9-34 

8.43 

.36 

6.76 

•19 

3.46 

•34 

6.47 

8 

3«T 

.13 

.38 

1.22 

2.93 

.38 

5.0T. 

3.18 

6.78 

9.88 

8.60 

•57 

6.88 

.*0 

4,70 

.54 

6.W 

9 

15.0 

.10 

.12 

1.12 

3.73 

• 38 

5.00 

3.18 

8.83 

9-69 

8. *8 

.56 

6.78 

.*0 

5.36 

•56 

6.71 

10 

3.7 

.37 

.38 

1.33 

2-93 

1.09 . 

6.68 

3.18 

14.13 

15.10 

12.59 

■6V 

10.07 

1.1V 

6.38 

.77 

9.25 

11 

3-3 

.27 

.52 

1.30 

2.79 

.76 

3.84 

1.77 

V.62 

5.49 

*•39 

.54 

5.2V 

.80 

V.71 

.49 

4.76 

12 

*.3 

•32 

.1* 

1.30 

3.07 

.96 

6.32 

1.77 

4.83 

5.63 

4.34 

-53 

5.18 

1.05 

3.66 

.50 

V.90 

13 

3-7 

.39 

-38 

1.30 

3-30 

1.29 

6.8V 

1.77 

V^7 

3. as 

V.46 

•52 

3.09 

1.3V 

6.72 

.51 

4.92 

1* 

3.7 

•39 

.59 

1.33 

2.92 

1.13 

6.74 

2.06 

8.30 

9.23 

T.49 

.62 

7.10 

1.48 

6.51 

.75 

6.51 

15 

9-5 


.00 

1.06 

3.58 

.06 

V.1V 



-- 


-- 



V.ll 



16 

6.2 


.09 

1.06 

3.35 

.01* 

3.89 


-- 

-- 



-- 


V.05 



17 

6.3 

•5S 

.*o 

1.33 

3.37 

1.90 

7.73 

1.68 

V.ll 

4.63, 

3.61 

.41 

3.13 

2.0V 

T.93 

.37 

4.60 

18 

1-9 

.46 

.1*3 

1.33 

3-lfl 

1.47 

7.13 

1.88 

V.10 

4.53 

3.69 

.V2 

5.27 

1.33 

6.97 

.37 

4.68 

19 

3.3 

.A 

.50 

1.28 

2.8V 

.68 

5.73 

.76 

e.w 

9 . 1 T 

7.7V 

.8V , 

11.72 

.68 

3-48 

.80 

11.1V 

20 

6-9 

.10 

.20 

I.16 

3.V2 

•3V 

4.96 

.76 

8.34 

9.17 

8.10 

.88 

12.27 

•3V 

V.8Q 

•8V 

U.69 
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s* 








BieoreticaJ. 


Center of pressure 


-50 
.29 
.27 
.*1 1.27 

.32 1.27 

.50 l.U 
.52 l.U 
.18 l.U 
.3* 1.27 

.33 1.27 

.26 1.21 
• 32 1.2T 

-31 1.27 

.61 l.ll 
-39 1.22 

-33 1-22 

.26 1.19 

.23 1.19 

.57 l.ll 
-Vo 1.31 
A3 1.31 
.23 1,19 

.er i.n 
.37 l.ll 
.63 1.11 


■25 1.19 
.28 1^9 
.33 1.25 
•37 1.25 
.33 1.22 
.30 1.22 
.39 1.22 
• 39 1.22 
.21 1.18 
.33 1.19 
.12 1.22 
•52 1.33 



>E<W) 

ic 

!.01 

2.U 

1.72 

!-oe 

2.66 

1.6*. 

.89 

8.51 

7.3* 

■92 

8.78 

7.11 




.61 
.62 
.62 
.58 
.11 
.11 
.11 
.60 
.60 
.62 
-58 
-57 
-59 
.13 
.67 
.13 
.11 
.71 

9-15 7.16 .79 

8.96 7.51 -79 

8.72 T-95 .66 

9.05 6.90 .76 

9.19 6.73 -71 

5-82 1.25 .17 

9.65 7.25 .76 

8.58 7.21 .76 

15.61 13.78 .81 

9.33 8.82 .86 

5-21 1.27 .26 

3.01 2.69 .26 

9.00 8.38 .56 

9.31 8.17 .56 

8.23 9-51 8.20 .55 

8.81 9.92 8.65 .58 

17.96 20.27 16.87 .70 

20.33 17.06 .71 

20.11 17.06 .71 

1.86 3.90 .15 

1.30 I 1.71 3-92 .15 


2.06 1.00 5.15 
2.06 1.01 5.62 
1.27 10.65 11.50 
2.06 15.97 17.25 
2.06 15.85 17.58 
2.0 6 15.97 17-53 
2.06 15.99 17.75 
2.06 9.38 12.02 
2.51 3-30 3.71 
2.31 8.29 8.95 

2-51 3-31 3-51 

2.51 3-36 5.82 

1.50 7.85 8.11 

1.09 8.51 9.15 

1.09 8.50 8.96 
8.30 8.72 


lift 


3(Cl tt ) s 


1.31 11.05 

7.97 16.29 


6 .87 15.12 
5-68 12.21 

7.88 16.93 
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Siender-body theory 

Strip theory 

Wing alone 
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Body -radius, wing -semispan ratio , r/s m 
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Figure 5 -Comparison of theoretical values of (Jr-) for 

1 c r 'W(B) 

triangular wing with no trading-edge sweep. 
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O 2 4 6 

Effective wing aspect ratio , /SA 
(a) No traiiing-edge sweep * 


Figure 6. - Center of pressure for wings according to linear 
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O 2 4 6 8 10 


Effective wing aspect ratio, /3 a 
( b) No midchord sweep . 

Figure 6.- Continued. 
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O .2 . 4 .6 .8 _ 1.0 

Theoretical center of pressure , (j)^ 


Figure 8. - Correlation between experimental and estimated 
centers of pressure for triangular wing-body combinations. 








Experimental center of pressure, (yh 


6M 
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Line of perfect, 
agreement 


□ afterbody 
no afterbody 


Theoretical center of pressure, (y)^ 


Figure 9.- Correlation between experimental and estimated 
centers of pressure for rectangular wing-body 
combinations. 
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o .2 4 .6 .1 8 1.0 

• t 

Theoretical center of pressure , ( j)q 

t 

Figure 10. - Correlation between experimental and estimated 
centers of pressure for trapezoidal wing -body 
combinations . 
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